Magnetoresistance of Three-Constituent Composites: Percolation Near a Critical Line 
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Scaling theory, duality symmetry, and numerical simulations of a random network model are used 
to study the magnetoresistance of a metal/insulator/perfect conductor composite with a disordered 
columnar microstructure. The phase diagram is found to have a critical line which separates regions 
of saturating and non-saturating magnetoresistance. The percolation problem which describes this 
line is a generalization of anisotropic percolation. We locate the percolation threshold and determine 
the values of the critical exponents tu — t± = sh = s± = 1.30 ± 0.02, v = 4/3 ± 0.02, which are 
the same as in two-constituent 2D isotropic percolation. We also determine the exponents which 
characterize the critical dependence on magnetic field, and confirm numerically that v is independent 
of anisotropy. We propose and test a complete scaling description of the magnetoresistance in the 
vicinity of the critical line. 



I. INTRODUCTION 

Magnetotransport in composite conductors has at- 
tracted increased attention due to the discovery of new 
effects, like the appearance of non-saturating magnetore- 
sistance in a metal/insulator columnar composite (de- 
noted by M/I), induced by the Hall effect in the M (i.e., 
the metallic) constituent n This means that, even when 
the M constituent has no intrinsic magnetoresistance but 
only a Hall resistivity, an induced magnetoresistance ap- 
pears in the mixture and continues to increase as B 2 
(B is the applied magnetic field) indefinitely whenever 
the Hall-to-Ohmic resistivity ratio of that constituent is 
much greater than 1. By contrast, in a columnar compos- 
ite of normal conductor and perfect conductor (denoted 
by M/ S), while there also appears an induced magnetore- 
sistance, it saturates when B is comparably large.EJ Re- 
lated new effects were also found in periodic microstruc- 
tures of either the M/I type or the M/S type, where 
the induced_magnetoresistance often exhibits a strong 
anisotropy.lm 

In this paper, we study the magnetoresistance of three- 
constituent composites with a random columnar mi- 
crostructure. Specifically, we analyze a composite con- 
sisting of cylindrical / and S inclusions, though not nec- 
essarily circular-cylinders, in an M host film, with cylin- 
der axes perpendicular to the film, and with both a mag- 
netic field and a uniform macroscopic or volume averaged 
electric current applied in the film plane. We assume 
that the M constituent has a Hall effect, with a Hall-to- 
transverse- Ohmic resistivity ratio H = pnaii/ Pohmic = 
/i|B | , where p is the Hall mobility. We are concerned with 
the effective resistivity tensor p e of this system at large 
H (p e is defined by (E) — p e (J), where (E) and (J) arc 
the volume averaged electric field and current density.u) 

This system might appear to be inherently three- 
dimensional (3D), because the Hall effect will generate 
an electric field with a component perpendicular to the 
film even with an in-plane applied current. However, the 



electric field perpendicular to the film plane vanishes be- 
cause of the columnar S. inclusions. Moreover, this 3D 
problem can be reducedtl to that of calculating the effec- 
tive conductivity of a two-dimensional (2D) composite of 
perfect insulator / with o~i = 0, perfect conductor (or su- 
perconductor) S with 175 = oo, and anisotropic 2D metal 
M with conductivity tensor 
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The conductivities o~m,\\ and <7m,_l correspond to the in- 
plane directions parallel and perpendicular to the applied 
magnetic field B. This transformation underlies our fur- 
ther discussion. The form assumed for &m means that 
the M constituent remains an isotropic conductor, even 
in a magnetic field. This excludes "open orbit" conduc- 
tors, but includes the possibility of some intrinsic mag- 
netoresistance: In that case both the transverse and the 
longitudinal Ohmic resistivities, pM and XpM, would de- 
pend upon B. Our subsequent discussion will assume, 
for simplicity, that A = 1 and pm are both independent 
of B, as would be the case if M were a free electron or 
free hole conductor. In that case, H is simply propor- 
tional to |B|. The three-constituent problem in 2D, but 
with a scalar &m, was treated previously in Ref. ||. Scal- 
ing and simulational studies, of magnetotransport in 3D, 
two-constituent M/I composites with a random inotropic 
microstructure, were also performed previously.^ 



II. MAGNETORESISTANCE AND 
PERCOLATION 

The relation of this problem to percolation can be 
understood by considering very large H. In this limit, 
o~m,-L 0, and within the metallic constituent the in- 
plane current prefers to flow in the o~m,\\ direction. For 
large enough volume fraction p$ of the S constituent, 
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and small enough volume fraction pj of /, there always 
exist current paths between S grains that are everywhere 
parallel to the high-conductivity direction in M [cf. Fig. 
0(a)]. Therefore, in this regime, at large H the effective 
composite conductivity (and hence resistivity) saturates 
at some finite value. In the opposite case (small ps, large 
pi), the current will be forced in some regions to flow in 
the low conductivity direction [cf. Fig. 0(b)]. In this case, 
the macroscopic or effective conductivity of the system 
will be proportional to o m,x! hence, the effective resis- 
tivity p e satisfies p e ~ 1/ctm,± ~ H 2 , i.e., p e will never 
saturate. 

Next, we qualitatively discuss the expected "phase di- 
agram" of this composite. The effective conductivity a e 
of the 2D mixture is a 2 x 2 tensor (like the effective 2D 
resistivity p e — cr" 1 ), whose components depend on the 
applied magnetic field. If the in-plane microstructure is 
isotropic, then the principal axes of a e are determined by 
B. At large H, both a e » and a e ± can either decrease 
as H~ 2 for arbitrarily large H, or else saturate at some 
finite value. In principle, there could be a region in the 
phase diagram where the resistivity saturates along one 
principal axis, but not along the other. But we expect 
no such region in an infinite system. This conclusion, as 
well as the location of the critical line, follow from the 
relation between the present problem and the anisotropic 
percolation model in 2D, as we now explain. 




FIG. 1. Schematic of transport in a 2D composite of a 
perfect conductor S, insulator /, and a very anisotropic metal 
M (the host medium in this figure). The current prefers to 
flow between S grains only in the direction of high metal- 
lic conductivity a M «, i.e., parallel to B. At high S volume 
fraction ps, such a path always exists, as shown schemati- 
cally in (a). At smaller ps or larger pi, the current must 
sometimes flow through the metal in the low-conductivity di- 
rection (<jm,± ~ H~ 2 ), as shown in (b). This behavior leads 
to a non-saturating effective resistivity. 

Anisotropic percolation is usually defined in terms of 
a random resistor network (RRN) , where the bond occu- 
pation probability depends on the bond orientation. In 
our 2D network model those probabilities will be chosen 
as pn = pm +ps for bonds along B, the high conductivity 
direction of the M constituent, and p± = ps for bonds 
perpendicular to B, the low conductivity direction. For 



an infinite network, the percolation threshold for both 
directions occurs whenQ 

1 = P\\ + Px = 2ps + Pm =>Pi = Ps- (2) 

This threshold separates the regimes of saturating and 
non-saturating magnetoresistance. It agrees with a pre- 
diction phased on the effective medium approximation 
(EM A), Em as well as with numerical results below. In a 
composite with less symmetry between the I and S con- 
stituents, the threshold will usually differ from pi = ps- 
The predicted phase diagram is shown in Fig. 0. For 
ps > 0.5, the S constituent percolates, hence the com- 
posite is a perfect conductor. Similarly, if pi > 0.5 then 
the combination of M and S constituents does not per- 
colate, hence the composite is a perfect insulator. These 
behaviors are independent of how much of the M con- 
stituent is present. The remainder of the phase diagram 
is divided into saturating and non-saturating regimes by 
the line pi = ps- In both regimes, our numerical results 
indicate that it is irrelevant whether the M constituent 
percolates by itself (as in the area below the dotted line), 
or whether only the aggregate of M and S constituents 
percolates. 
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FIG. 2. Schematic phase diagram for an M/I/S com- 
posite with random columnar microstructure. pi and ps are 
the volume fractions of / and S; the M volume fraction is 
Pm — 1 — pi —ps- In the shaded regions, either the S or the / 
constituent percolates. The critical line pi = ps separates re- 
gions with saturating and non-saturating magnetoresistance. 
The dotted line corresponds to percolation of the metallic 
constituent alone (pm > 0.5 or pi + ps < 0.5); no critical 
behavior appears to be associated with this line. 



III. CRITICAL EXPONENTS AND SCALING 
PROPERTIES 

The critical behavior on approaching the critical line is 
described by a number of critical exponents. One of these 
governs the two correlation lengths, £|| and which 
diverge at the critical line. But we expect that the di- 
vergence of both will be governed by the same critical 
exponent v, because their ratio £ii/£± is determined by 
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remains finite on crossing the critical line 
write 



iff 



Thus we 



(3) 



Moreover, since the correlation length is just a geomet- 
rical property of a percolating system, it should have 
the same value v a as in anisotropic percolation, i.e., 
v — v a . An earlier renormalization group analysis pre- 
dicts that anisotropic, percolation is governed by the 
isotropic fixed point .113 We therefore expect that v a has 
the same value as in 2D isotropic percolation, wiuch 
has been shown both analytically^ and numericallytJ to 
equal 4/3. However, a similar proof has not been given 
for anisotropic percolation;ll3 moreover, some time_ago it 
was reported that v does depend on anisotropyJl3 Our 
numerical results below give v = 4/3 ±0.02, independent 
of anisotropy. 

The magnetoresistance is described by other critical 
indices. At any point such that pi < ps (the region 
below the critical line in Fig.|^), the effective conductivity 
saturates at some finite value (i = || or i =_L): 



g sat ,i = lim a sa t,i{H). 



(4) 



As we approach the critical line (say, along the line pu — 
const), g sa t t i must tend to 0, since on the other side of 
that line the conductivity vanishes at large H as H~ 2 . 
Hence we can introduce the critical exponents tu and t± 
according to 



9sat,i ~ (PS - PlY 



(5) 



In the region pi > ps, the effective conductivity is pro- 
portional to l/H 2 . Therefore, we can define the finite 
limiting value 



Qnon — sat,i — lim H (7 non - 
H—>oc 



sat,i{H). 



(6) 



Since g n on-sat.i must diverge as the critical line is ap- 
proached, we can define the critical exponents sm and 
sj_, according to 



9n 



-sat,i ~ (Pi —Ps)~ 



(7) 



Just as in the case of is, we expect sii = sj_, tu = t±. We 
might expect that the actual values of t and s are deter- 
mined, once again, by the connection to anisotropic per- 
colation. However, this connection is more subtle than 
for v, because t and s refer to different physical quanti- 
ties in the present problem, than in anisotropic percola- 
tion. Nonetheless, it is reasonable to expect that both 
problems belong to the same universality class, implying 
Si = ti = 1.30, the values of the two conductua,ti£-£xpo- 
nents in conventional 2D percolation problems.EZHijNote 
also that EMA yields Sj = = ljffl and that continuum 
composites sometimes belong to a different universality 



class of percolation, with different values of the critical 
exponents, depending upon the microstructural details.EHI 
Finally, consider systems exactly at the percolation 
threshold, i.e., at pi — ps- Any finite-size sample 
inevitably falls into either the percolating class, with 
p e ~ H° as H — > oo, or the non-percolating class, with 
p e ~ H 2 . On increasing the size of the system, one finds 
that the field where a saturating sample achieves satu- 
ration also increases. The same size effect also holds for 
non-saturating samples: the characteristic field at which 
p e starts to vary as H 2 increases with the size of the 
system. Below that field, or at any field in the limit of 
a system of infinite size, the behavior of the magnetore- 
sistance is described by yet another critical exponent <$j, 
defined by 



\H[ 



(8) 



The EMA analysisBfl yields 6, t = l. 

In fact, the result S j = 1 also follows from a duality 
argument, if one assumes that 5\\ = S±. The duality 
principles gives 
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/ 1 



1 



1 1 



\0M_L a M\\ a S 07 



It immediately follows that the anisotropic percolation 
thresholds for the two directions ||, J_ must coincide, for 
otherwise the singular behaviors of <j e n and a e< ±_ would 
not cancel, as required by this equation. It also fol- 
lows, rigorously, that = s±_ and t±_ = sii. Using 
the homogeneity of c e i|, cr ej ± as functions of the vari- 
ous constituent conductivities, and noting that <tj = 
and a$ — oo, we can rewrite this equation as 
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Note that the relation pi — ps for the percolation 
threshold follows rigorously from this equation if the mi- 
crostructure is invariant under the interchange of the S 
and / constituents. At the percolation threshold, and 
when \H\ ^> 1. this reduces to 
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The assumption that Su = S±, and hence the final re- 
sult (5|| = S± = 1, are supported by the physical picture 
of the anisotropic percolation process, and is consistent 
with the expectation that t|| = tj_ and S|| = sj_. These 
equalities also follow if a simple scaling description ap- 
plies to both a e> \\ and cr ej j_ with the same scaling variable. 
Thus, for an M/I/S columnar composite precisely at the 
critical composition (i.e., on the phase boundary line be- 
tween saturating and non-saturating magnetoresistance), 
we expect to find 
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Near the transition, where \pi — Ps\ ^ 1 an< i 
cm,_l/cm,|| — l/-ff 2 *C 1, we expect that a scaling de- 
scription of the critical behavior is applicable. In view of 
the preceding discussion, that description can be formu- 
lated as follows: 

^-signAplApr'F^Z); (9) 

Pm,\\ 

A P = Pl -p Sl i=\\,±, Z= |F||Ap|*signAp. 

The scaling functions Fi(Z) should have the following 
asymptotic forms for extreme values of the scaling vari- 
able Z: 

!-Z° for Z < -1, i.e. p s > p u 
Z 2 for Z > 1, i.e. p s < pi, (10) 
Z for \Z\ < 1, i.e. p s < pi- 

The first two lines in this expression follow from the fact 
that we expect to have p e .i ~ -ff in the saturating regime 
and /? e ,i ~ H in the non-saturating regime. The third 
line results from the necessity to cancel the dependence 
of p e ^i upon Ap when Ap — > 0. This kindrpf behavior 
was already found previously using EM A, Bel where the 
incorrect value t = 1 was found, as is usual when that 
approximation is invoked. Note that, even though Fi(Z) 
has a different analytic form for large Z 7 depending on 
the sign of Z or Ap, these functions are expected to vary 
smoothly when Z passes through 0. 

Finally, note that duality is a general symmetry prop- 
erty of 2D systems, and does not require that the I and 
S inclusions have similar shapes or spatial distributions. 
Thus, the conclusions regarding values of s, t, S are valid 
even if the line of critical compositions differs from the 
simple straight line ps — pi, which is valid only for the 
case where the I and S constituents appear in the com- 
posite in a symmetric fashion. 

IV. RESULTS 

We carried out calculations on simple-square-lattice 
RRN's, choosing the resistors in accordance with the 
model described above. The calculation,-was done us- 
ing the Y-A bond elimination algorithmic Besides be- 
ing very efficient in 2D, this algorithm allows the inclu- 
sion of both perfectly insulating and perfectly conduct- 
ing bonds without any approximations. This technique 
is much more efficient than the techniques available for 
3D networks. That is why the results obtained here for 
columnar systems are much more detailed and accurate 
than the results obtained previously for 3D isotropic-sys- 
tems by simulations of 3D random network models 

The networks for our results were constructed as fol- 
lows: each bond was independently and randomly chosen 
to be insulating, perfectly conducting, or metallic, with 
appropriate probabilities. A metallic bond was assigned 
an appropriate conductance, depending on whether it 



was oriented parallel or perpendicular to the magnetic 
field. To calculate g sa t,i, as can be seen from Eqs. (|l|) 
and (||) , the conductances of the constituents should be 
taken as 

07 = , a s = oo , a M ,\\ = 1 , ctm.jl = 0. (11) 

But in order to calculate g n on-sat,i, one needs to multiply 
all the conductances by H 2 and take the limit H — > oo, 
which leads to 

07 = , a s = oo , cr M .|| = oo , a M ,± = 1. (12) 

The definitions (||) and (fy for the exponents s and t are 
exactly valid only for an infinite system. Instead of using 
these definitions directly, we adopted the finite-size scal- 
ing approach, as described below. For a finite-size system 
generated exactly at the percolation threshold pi = ps, 
the system size L x L is always smaller than the (infi- 
nite) correlation length £. From Eq. (||) it follows that 
such a system behaves like a system at another volume 
fraction satisfying \pj — ps\ ~ L -1 '". Therefore, the 
average of g sa t,i or g n on-sat,i over many such systems 
should scale as (g sa t,i) ~ L~ u/u or (g n on-aat,i) ~ L s ^ u . 
(Due to the appearance of percolating samples, the lat- 
ter average is actually infinite. A finite result is obtained 
by averaging only over non-percolating samples. This 
procedure does not change the scaling form, because at 
the percolation threshold the fraction of samples which 
are non-percolating is asymptotically size- independent, 
as discussed below). 

Using RRN's with L ranging from 100 to 2000 and 
Pi = Ps, we estimated Sj and U for pu = 0.2, 0.5, 0.8, 
and, with less accuracy, for several other values of pu- 
(At pm = 0.8 we could estimate only t\\ and sj_, since 
these systems usually percolated only in one direction but 
not in the other.) The finite size scaling assumption is 
supported very well, as demonstrated in Fig. |^ for pm — 
0.2. Assuming v = 4/3, we find for all values of pm 
studied that tu = t± = sn = sj_ = 1.30±0.02. This value 
supports the hypothesis that the problem belongs to the 
same universality class as both isotropic and anisotropic 
percolation in 2D. 

The assumption that v = 4/3 can be tested for our 
problem by considering percolation in finite-size systems. 
For such systems, the "percolation threshold" in the di- 
rection i is naturally defined as the value of ps for which 
exactly one half of the systems percolate in that direction: 
Pi = N sa t,i/N tota i = 1/2, where P, is the probability of 
finding saturating behavior in the direction i. Fig. ^ then 
shows that the percolation threshold depends on the sys- 
tem size L. In fact, a finite system has two "percolation 
thresholds", psq.± andpso,||; located respectively above 
and below the percolation threshold of an infinite system, 
and apnroaching pi as \pso.i ~Pi\ ~ L^ 1 ^ with increas- 
ing Lr3 We have checked these hypotheses numerically 
by calculating \pso _l — PSo || I f° r f° ur linear sizes L at 
p M = 0.2 and 0.5^ and verified that v = 4/3 ± 0.02, 
independent of anisotropy. 
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FIG. 3. Calculated (g sa t,i) and (g no n-sat,i) [cf. definitions 
(^) and (|^)] plotted on a log-log plot versus linear size L of 
the RRN for pu = 0.2 and pi — ps — 0.4. Error bars are 
smaller than the dot sizes. The slopes of the two lines on this 
plot equal ti/v and —Si/v, and are consistent with the value 
0.975, corresponding to t = s = 1.30 and v = 4/3. 

Fig. [| also makes clear that the fraction of samples 
which are non-percolating is size-independent at the per- 
colation threshold. Specifically, Fig. |] shows the fraction 
of percolating samples of a given size vs. concentration 
ps. The lines corresponding to different sample sizes L 
intersect at the percolation threshold pj = p$ (vertical 
line). Thus, as L — > oo, the fraction of the samples perco- 
lating at pj — ps approaches a certain limiting value de- 
termined by the other parameters of the system, namely, 
the direction of percolation and the concentration pu of 
the metal. This limiting value can be related to the span- 
ning probability in the case of isotropic percolation in a 
system of rectangular shape, which swas found to depend 
on the aspect ratio of the rectangle.Ej In the present pa- 
per, we only consider systems with the aspect ratio equal 
to unity; however, the correlation lengths are different in 
the two directions. Thus, may be thought of as an 

effective aspect ratio for the present problem. 

In order to test the scaling behavior at finite values 
of H and Ap, we also simulated random networks at 
Ps = Pi with large but finite values of both H and L. 
Since the RRN's required for such calculations contain 
bonds with widely different conductivities [cf. Eq. ([[])], 
special care must be taken when performing these calcu- 
lations on large systems. c3 For that reason we wrote spe- 
cial code for storing numbers as large (small) as roughly 
10 ±6x10 . Using this code, we studied RRN's of size up 
to 4000 x 4000 in fields ranging from up to H 2 = I0 8 . 
One useful way to exhibit those results is as plots of 
ln[(H 2 + l)o- e ,i] vs. ln(H 2 + I) for different values of the 
linear size L — see Figs. || and ||. Fig. ^ shows <r e m and 
<7 e j_ vs. magnetic field in systems with sizes ranging from 
100 to 4000. For the larger systems, one can clearly see 
a "critical" range of fields, in which the magnetoresis- 
tance is consistent with the scaling form of Eq. (||) with 
6±_ = <5|| = 1.00 ± 0.07. For fields within and below that 
critical range, the behavior of o~ e ^ is independent of L, 
and is also independent of whether its value does or does 



not saturate at higher fields, as demonstrated in Fig. [|. 
However, the intercepts of those linear dependencies are 
different for the different directions: o- e; i ~ (\H\ — H ,i) 5i , 
with positive i?o,_L and negative -Ho,||- This is due to the 
fact that the system always conducts better in the di- 
rection parallel to the applied magnetic field, as can be 
surmised from Eq. (|l|). Figs. || and || bear out the forms 
hypothesized earlier for Fi(Z) — see Eq. (|l0|). In particu- 
lar, Fig. U bears out the expectation that, for small \Z\, 
Fi(—Z) = —Fi(Z), and thus that the scaling functions 
are smooth at Z = 0. Figs. |5| and ^ also show that both 
cr e ii and <Te,J- are independent of L when \H\ <C L l l v ', in 
agreement with Eqs. (Eol) and (O) below. 
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FIG. 4. Fraction Pi = N sa t / Ntotai of samples which sat- 
urate in the direction i, plotted versus 5 volume fraction ps, 
at fixed M volume fraction pu =0.2, for two different sizes 
L. The vertical line denotes the percolation threshold in an 
infinite system, as predicted by Eq. (^). For a finite system, 
the "percolation threshold" is the point where one half of the 
samples percolate: Pi = 0.5. Such thresholds are unequal in 
the directions parallel and perpendicular to B, but approach 
the same infinite-size value, ps = 0.4, as L — > oo. 

In order to find the form of the scaling functions Fi(Z) 
[defined by Eq. (^] from these numerical results, it is coft. 
venient to also invoke the finite-size-scaling hypothesis:EJ 
In a system of finite linear size L, when the correlation 
length £ satisfies £ 3> L, \Ap\ should be replaced by 
Ci/L 1 ^ in all the expressions of Eq. (||): 



Pe.i 

Pm. ii 



£*/"signApFi(Z), Z = signAp\H\/ . (13) 



Note that the constant C\ has been absorbed into the 
definitions of Z and Fi(Z). The sign of Ap, appearing in 
Eq. (H|), should now be determined from the actual be- 
havior (i.e., percolating vs. non-percolating, or saturating 
vs. non-saturating, in the direction i) of each particular 
sample. 

The plots of L- t l v p^i/ p M ,\\ vs. \H\L- l l v , shown in 
Fig. 0, clearly demonstrate that the results obtained for 
Pe.i, using different values of H and L, collapse onto a 
single curve (for a given i) when scaled in accordance 



5 



with Eq. flU). These fi gures constitute a quantitative 
graphical representation of the scaling functions Fn (Z) 
and F±(Z) for the special case pu = 0.5. Note that these 
two scaling functions appear to have a similar shape, up 
to a constant coefficient. That is qualitatively consistent 
with the results of EMA, which lead tcfl 
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1 + Pm 
1 - Pm 



In the case of the systems featured in Fig. [?], where pu — 
0.5, this ratio should be 9 according to EMA. By contrast, 
the simulation results plotted in that figure indicate that 
this ratio is between 4.5 and 5. As stated earlier, such 
quantitative discrepancies should come as no surprise, in 
view of the known deficiencies of EMA in the critical 
region near a percolation threshold. 
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govern the dependence on magnetic field in the critical 
region. The existence of scaling behavior in that region 
was confirmed and the forms of the scaling functions were 
obtained numerically. 
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FIG. 6. ln[(tf 2 + l)a e ,i] vs. \n(H 2 + 1) for systems of 
size L — 100 (triangles), 1000 (squares) and 4000 (circles), 
with pm = 0.5, pi — ps = 0.25. Open symbols correspond to 
a Ei x and the filled symbols to a e ii. The lines are drawn as 
guides for the eye. At large fields, all systems exhibit either 
saturating (H 2 ai ~ H 2 ) or non-saturating (H 2 <Ji ~ 1) mag- 
netoresistance. In the "critical" range of fields, clearly seen 
for the larger systems, the magnetoresistive behavior is given 
by Eq. (§) with 5± = 5y = 1.00 ±0.07. Note that, for the sake 
of visual clarity, we include only those systems in which a ei ± 
does not saturate and a e ^ does saturate [cf. Fig. |H). 



FIG. 5. ln[(# 2 + l)a e ,i] vs. \n(H 2 + 1) for systems with 
L — 100 and pm = 0.5 at the percolation threshold pj = ps- 
Averaging over different realizations is performed separately 
for saturating and non-saturating samples. Open symbols 
correspond to a sat ,± and cr n on.-sat,x and the filled symbols to 
°sat, II and CTnon-sat, ||- F° r values of H within and below the 
"critical" range of fields (as defined in the text), the conduc- 
tivity in a given direction does not depend on whether or not 
it saturates at stronger fields H. 

V. DISCUSSION 

We have demonstrated that the phase diagram of an 
M/I/S random columnar composite in an in-plane mag- 
netic field exhibits a critical line between saturating and 
non-saturating regimes of magnetoresistance. The be- 
havior near this critical line was interpreted using an 
anisotropic percolation model. The critical exponents 
t = s and v were found, numerically, to be the same 
as in the more conventional 2D isotropic random perco- 
lating networks. Also found were the exponents which 



The critical line discussed here could be studied experi- 
mentally using a doped semiconductor film as the M con- 
stituent, with a random collection of etched perpendicu- 
lar holes as the / constituent, and a random collection of 
perpendicular columnar inclusions, made of a high con- 
ductivity normal metal, playing the role of S. Extremely 
low temperatures or very clean single crystals would not 
be required in order to observe this critical line. What 
would be necessary is a large contrast at each stage of the 
following chain of inequalities ps <C pu "C H 2 pM <§C pi- 
ll Si-doped GaAs is used as the M host, with a negative 
charge carrier density of 1.6 x 10 18 cm -3 and a mobility 
p = 2500 cm 2 /Vs at a temperature of 90 K, as in the 
experiment described in Ref. |^, then a magnetic field of 
40 Tesla would result in H = —10. Such a material would 
have an Ohmic resistivity of 1.6 x 10~ 3 f2cm, about 1000 
times greater than that of Cu. Thus, using Cu for the S 
inclusions and etched holes for the I inclusions, all the 
above inequalities could be satisfied without difficulty. 

This is apparently the first experimentally acces- 
sible and technologically promising system in which 
anisotropic percolation is relevant. Our numerical results 
for the exponents v, t, s, and 5 are consistent with the 
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assumption that this problem belongs to the same univer- 
sality class as the usual 2D isotropic percolation problem, 
and confirm that v is independent of anisotropy. 
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FIG. 7. Plots of Fi(Z) = sigaApL- t/l, pe ti /p UA] vs. 
Z = signAp \H\L~ l l u for networks with pm = 0.5 and vari- 
ous large but finite values of H and L. The top frame shows 
Fi(Z) in the non-saturating regime, when Z > 0, while the 
bottom frame shows Fi(Z) in the saturating regime, when 
Z < 0. The shapes of the points correspond to the value of L, 
as in Fig. ^. Open symbols correspond to F±(Z), while filled 
symbols correspond to F\\{Z). The horizontal lines surround- 
ing the symbols are parts of error bars. Inset: expanded view 
of Fi(Z) for small values of \Z\ and both signs of Z, showing 
the smooth linear behavior of Fi(Z). 
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